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The agent responsible for the recent severe acute respiratory syndrome 
(SARS) outbreak is a previously unidentified coronavirus. While there is 
a wealth of epidemiological studies, little if any molecular characteri- 
zation of SARS coronavirus (SCoV) proteins has been carried out. Here 
we describe the molecular characterization of SCoV E protein, a critical 
component of the virus responsible for virion envelope morphogenesis. 
We conclusively show that SCoV E protein contains an unusually short, 
palindromic transmembrane helical hairpin around a previously unidenti- 
fied pseudo-center of symmetry, a structural feature which seems to be 
unique to SCoV. The hairpin deforms lipid bilayers by way of increasing 
their curvature, providing for the first time a molecular explanation of 
E protein’s pivotal role in viral budding. The molecular understanding 
of this critical component of SCoV may represent the beginning of a con- 
certed effort aimed at inhibiting its function, and consequently, viral 
infectivity. 

© 2004 Elsevier Ltd. All rights reserved. 
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Introduction 


A previously unidentified member of the corona- 
viridae is the etiologic agent responsible for the 
recent severe acute respiratory syndrome (SARS) 
outbreak.’ Viral genome sequencing, combined 
with protein phylogenetic analyses, have shown 
that SARS coronavirus (SCoV) belongs to a new 
subfamily within the coronaviridae.”* In an effort 
to better understand the components contributing 
to the pathogenic mechanism of the virus, we 
have structurally analyzed the transmembrane 
domain of SCoV E protein using Fourier Transform 
Infrared (FTIR) spectroscopy, X-ray scattering, 


Abbreviations used: SARS, severe acute respiratory 
syndrome; SCoV, SARS coronavirus; MHV, mouse 
hepatitis virus; TGEV, transmissible gastroenteritis virus; 
BCoV, bovine coronavirus; IBV, infectious bronchitis 
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ATR, attenuated total reflection; FTIR, Fourier Transform 
Infrared. 
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electron microscopy and global searching molecu- 
lar dynamics simulations. 

Functionally, coronavirus E protein’s pivotal role 
in viral morphogenesis has been studied exten- 
sively. Co-expression of Mouse Hepatitis Virus 
(MHV) E and M proteins on their own has been 
shown to result in the production of virus-like 
particles.*° This fundamental observation proved 
that neither the nucleocapsid nor the viral spike 
are needed for viral budding. Work on other 
viruses, such as Transmissible Gastroenteritis 
Virus (TGEV),° Bovine coronavirus (BCoV)*® and 
Infectious Bronchitis Virus (IBV)’ has shown this 
to be a general phenomena of the coronaviridea. 
The interaction between the cytoplasmic domains 
of the two proteins is thought to take place in 
pre-Golgi compartments.* 

The importance of the role played by E protein 
in virus budding and morphogenesis is further 
strengthened by the observation that in a number 
of coronaviruses, expression of M protein on its 
own does not produce virus-like particles.7~”” In 
contrast, expression of IBV or MHV E protein on 
its own causes the release of vesicles containing 
E protein, thereby pointing to the importance of 
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Figure 1. Top panel: sequence of the SCoV E protein 
(strain TOR2 NP_828854). Sequence of the synthesized 
peptide is boxed, while gray shading indicates the extent 
of the hydrophobic stretch of the protein. The iodinated 
Phe residue is indicated by dark shading. Bottom panel: 
multiple sequence alignment of the hydrophobic region 
of E proteins belonging to the four different corona- 
viridea subgroups. Gray shading indicates the hydro- 
phobic amino acid residues according to the GES scale. 
Sequence codes are given in the Materials and Methods. 


coronavirus E protein in the budding process.”” 
Interestingly, cells infected with MHV release 
vesicles containing E protein as well.’ Finally, 
mutations in MHV E protein cause marked mor- 
phological changes in the resulting viruses.’ 


In the cell, evidence has been accumulating 
suggesting that IBV E protein localizes to the 
Golgi apparatus,’ the site from which corona- 
viruses are known to invaginate into mature 
virions. During the expression of E protein the 
Golgi apparatus changes its morphology 
dramatically," which could explain in part E 
protein’s ability to induce apoptosis.’*!% 

Structurally, E proteins are conserved within the 
different coronaviridea groups, yet exhibit little 
sequence similarity between the groups (Figure 1). 
In line with its belonging to a new group within 
the coronaviridea, SCoV E protein differs substan- 
tially from other coronavirus E proteins.* Whether 
this has any bearing upon the marked difference 
in pathogenicity of SCoV versus other corona- 
viruses remains to be seen. Finally, SCoV E protein 
shares no significant sequence similarity with any 
other known protein. 

It is possible, however, to make the following 
generalizations regarding the sequence of all coro- 
navirus E proteins: they are small proteins (~75 
residues), with an unusually long hydrophobic 
stretch (25-30 residues) located between hydro- 
philic N and C termini (~8 and ~40 residues, 
respectively). SCoV E protein shares the same over- 
all characteristics, although only 20% identical 
with other coronaviridea E proteins:’ 74 residues 
in length with a stretch of 26 hydrophobic amino 
acid residues (Figure 1). 

The length of the hydrophobic segment of SCoV 
E protein has posed a problem with respect to 
assigning a topology to the protein (Figure 2(a)). 
The average length of a transmembrane a-helix is 
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Figure 2. (a) Hydrophobicity analysis of the SCoV E protein according to the GES scale” using a window size of 15 
and a hydrophobicity cutoff of AGwateroil = —27 kcal/mol.” The region of the synthesized peptide is boxed, while 
gray shading indicates the extent of the hydrophobic stretch of the protein. The location of the iodinated Phe residue 
is indicated by dark shading. (b) Inversion diagram of the SCoV E protein hydrophobic region outlining identity and 
similarity between the two inverted sequence elements, in continuous and broken lines, respectively. Hydrophobic 
amino acid residues according to the GES scale” are indicated by gray shading. 
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21 residues,* far shorter than the hydrophobic 
stretch of E proteins. In contrast, if a hairpin were 
to be formed by the long hydrophobic stretch, 
the resulting helices would be much shorter than 
the average. 

Protease digestion of MHV does not affect E 
protein’s molecular mass, suggesting that no 
part of the protein protrudes beyond the viral 
membrane." Antibodies raised against the N 
terminus of the protein could detect the protein 
only after treatment of the cells with Triton X-110, 
but not with digitonin.’ Both of the above studies 
concluded that the protein traverses the lipid 
bilayer once, with its C terminus located in the 
viral interior.” More recent results with an epi- 
tope tagged MHV E protein are indicative of the 
protein traversing the lipid bilayer twice, whereby 
both termini of the protein reside in the virus 
lumen.!° 

In the current study we present the results of a 
detailed structural examination of SCoV E protein, 
in which we determine the topology of the protein 
and the effects upon the lipid bilayer thereof. The 


results conclusively show that SCoV E protein 
has a highly unusual topology, consisting of a 
very short transmembrane helical hairpin. The 
hairpin forms an inversion about a previously uni- 
dentified pseudo-center of symmetry which seems 
to be unique to SCoV. Electron microscopy studies 
indicate that the protein dramatically distorts lipid 
bilayers, causing tubulation. Taken together, the 
structure of the protein facilitates a molecular 
explanation of the effects upon the bilayer curva- 
ture underpinning the function of the protein 
in vivo and may represent the beginning of a con- 
certed effort aimed at inhibiting its function and 
consequently viral infectivity. 


Results and Discussion 


E protein is decidedly helical 


In order to examine the secondary structure of 
the transmembrane domain (TMD) of SCoV E pro- 
tein, peptides were made which encompassed the 
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Figure 3. (a) Transmission FTIR spectrum of the SCoV E protein TMD reconstituted in DMPC lipid bilayers. (b) ATR- 
FTIR spectra of SCoV E protein TMD reconstituted in lipid bilayers, focusing on the amide I and amide II regions. 
Spectra were obtained after flushing air saturated with H,O (continuous line) or D,O (broken line) over the sample. 
Both spectra were normalized so as to adjust the absorption of the amide I mode to 1.0 OD. (c) ATR-FTIR spectra of 
SCoV E protein TMD reconstituted in lipid bilayers, focusing on the amide I region. Spectra were obtained with paral- 
lel (continuous line) or perpendicular (broken line) polarized light. (d) ATR-FTIR spectra as in c but focusing on the 
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entire hydrophobic region of the protein (Glu7- 
Arg38), as depicted in Figure 1. The peptide was 
purified and reconstituted into dimyristoylphos- 
phocholine (DMPC) lipid vesicles, to be used 
for FTIR spectroscopy. The transmission FTIR 
spectrum of the resulting proteo-liposomes focus- 
ing on the amide I vibrational mode (1600- 
1700 cm™') is presented in Figure 3(a). Since the 
amide I vibrational mode arises mostly from the 
peptidic C = O stretch, its frequency is well corre- 
lated with the secondary structure of the protein.'® 
Peaks resonating at 1655cm™', 1630cm™' and 
1645 cm~', correspond to a-helical, B-strand and 
random coil peptide segments, respectively.'® The 
amide I mode of SCoV E protein TMD is centered 
at 1657 cm™' with a peak-width at half height of 
23cm ~', both indicative of the very high helical 
content of the protein, and the absence of any 
other secondary structure component.”® 


E protein contains 26 residues embedded in 
the lipid bilayer 


FTIR spectroscopy can be used to delineate the 
extent of membrane incorporation of  trans- 
membrane proteins. This is achieved by observing 
the reduction in any vibrational mode containing 
significant contributions from the amide proton 
(e.g II mode:peptidic N-H deformation) upon 
exchanging the solvent from H,O to D,O. The 
reason is that the lipid bilayer protects the peptide 
from exchange, whereas amide groups exposed 
to the solvent undergo H*/D* exchange. Any 
vibrational mode containing the “new” N-D 
group will resonate elsewhere and H*/D* 
exchange can be quantitated by measuring the 
reduction in the amide II mode directly.’® 

As shown in Figure 3(b), upon flushing the 
membrane with air saturated with D,O for several 
hours, the reduction in the amide II peak, centered 
at 1545 cm~', is minimal. Since the amide I peak at 
1657 cm ' is not expected to undergo any intensity 
changes upon exchange from H,O to D,O, it is 
possible to calculate the extent of exchange by 
normalizing both spectra (in H,O and D,O) on the 
amide I peak. The resulting exchange rate of 
19 + 2% is indicative of only six out of 32 residues 
in the peptide undergoing Ht/D* exchange. In 
other words, 26 residues of SCoV E protein are 
protected from exchange as a result of their being 
embedded in the lipid bilayer. Note that this value 
corresponds exactly to the length of the hydro- 
phobic stretch of the protein. 


The transmembrane helical segment of E 
protein is virtually perpendicular to the 
lipid bilayer 


After establishing the helicity of the peptide and 
determining that most of it is embedded in the 
lipid bilayer, experiments were undertaken to 
determine its orientation using polarized FTIR 
spectroscopy of oriented membranes. Attenuated 


total internal reflection (ATR-FTIR) — spectra 
obtained using both parallel and perpendicular 
polarized light are depicted in Figure 3(c). The 
dichroism (ratio of absorption of parallel and per- 
pendicular polarized light) for the amide I peak 
(centered at 1657cm~') is 3.76 + 0.52. From this 
dichroism it is possible to calculate an exception- 
ally high order parameter for the peptide of 
S=0.82 +£0.17, a value nearing the theoretical 
limit.'7 When taking into account sample disorder, 
it is possible to convert the maximal value of the 
measured order parameter to a maximal tilt angle. 
Thus, the tilt angle of SCoV E protein helical 
elements from the membrane normal must be less 
than or equal to 6.3°. 


SCoV transmembrane domain is a 
helical hairpin 


So far three lines of evidence were obtained 
characterizing the transmembrane domain of 
SCoV E protein: (i) it is highly helical, (ii) it is 
composed of 26 amino acid residues and (iii) the 
helical elements are oriented normal to the mem- 
brane plane. Thus, there are two ways in which to 
describe the structure of the protein: either as a 
single long transmembrane helix, or as two short 
helices forming a transmembrane helical hairpin. 
Both models however, incur a hydrophobicity mis- 
match between the protein and the bilayer. In the 
case of a single helix, the hydrophobic stretch is 
too long relative to the bilayer thickness.+ On the 
other hand, the helical hairpin would have to be 
comprised of two very short helices in order to 
enable the 26 residues to traverse the membrane 
twice. Since the smallest possible loop is of 
three residues, each helix cannot contain more 
than 11-12 hydrophobic residues. This value is 
much smaller than the average of 21 residues per 
transmembrane helix." 

In order to conclusively determine the topology 
of SCoV E protein we iodinated Phe23, which is 
located at the center of the hydrophobic stretch of 
the protein (see Figure 1). FTIR spectroscopy has 
shown that the iodinated and unlabeled peptides 
are indistinguishable in any of the experiments 
described above (data not shown). Electron density 
profiles obtained using X-ray scattering should 
be able to pinpoint the position of the labeled 
electron-dense iodine. As shown in Figure 4, the 
electron density of lipid vesicles containing iodin- 
ated and unlabeled SCoV E protein transmem- 
brane domain exhibit normal density profiles.’ 
However, upon subtracting the densities form one 
another, clear density resulting from the iodine 
label is identified, located 16.5 A from the bilayer 
core, close to the lipid head-group region. Thus, 


+ Note that a single long hydrophobic helix cannot 
compensate for its elongated hydrophobicity by titling, 
since we have shown that the helix has a maximum tilt 
angle with respect to the bilayer normal of less than 6.3°. 
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Figure 4. Electron density profile (on arbitrary scale around the mean density) of the lipid bilayers obtained at molar 
protein-to-lipid (DMPC/SCoV E protein) ratios of P/L = 0.1, for iodinated (filled symbols) and non-iodinated (open 
symbols) protein, as obtained from Fourier synthesis of measured integrated Bragg peak intensities. The difference 
between the two curves indicates the position of the one iodinated phenylalanine per protein located adjacent to the 
head-group region, close to the hydrophilic/hydrophobic interface. For comparison, the electron density profile of 
the Influenza A M2 H* channel transmembrane domain, a tetrameric helical bundle,’ is depicted in broken line. 


one can clearly state that the labeled Phe23 is 
located adjacent to the head-group region, despite 
the fact that it is in the middle of the hydrophobic 
segment of the protein. Taken together, our data 
prove conclusively that the transmembrane 
domain of SCoV E protein forms a hairpin with 
unusually short hydrophobic helices. 


SCoV E protein contains an unusual 
palindromic symmetry 


After establishing that SCoV E protein forms a 
helical hairpin, close inspection of its sequence 
reveals an unusual feature. As shown in Figure 
2(b), inverting the protein about itself at Phe23 
indicates that the protein is virtually symmetrical 
with respect to the inversion point, thereby form- 
ing a palindrome. The only deviation from palin- 
dromic symmetry is located at Asn15, the 
significance of which is discussed below. 

In order to estimate the statistical significance of 
forming a palindrome within the helical hairpin 
framework, a Monte-Carlo simulation was per- 
formed. The sequence of each helix forming the 
hairpin was randomized, followed by assessing 
the extent of palindromic symmetry. The results 


indicate that the degree of inverted identity 
between the two helices of SCoV E protein (thereby 
forming a palindrome) is a rare event, observed at 
a frequency of ca 10°°. This would explain why to 
our knowledge palindromic symmetry in helical 
hairpins has not been observed before. 


E protein dramatically effects lipid bilayer 
structure and morphology 


The SCoV E protein contains the shortest trans- 
membrane hydrophobic helical hairpin known to 
date. Its 26 hydrophobic amino acid residues not 
only traverse the lipid bilayer as two helices, but 
must form a connecting turn between the helices. 
An immediate question then arises: What affect 
does the unusual structure of SCoV E protein 
transmembrane domain have upon the lipid 
bilayer? In order to answer the above question we 
undertook two lines of experiments, aimed at 
(i) examining the molecular structure of the lipid 
in detail, alongside (ii) establishing the effects 
upon the membrane morphology as a whole. 

The order of the lipid bilayer acyl chains was 
examined using polarized FTIR spectroscopy of 
oriented multilayered DMPC vesicles containing 
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Figure 5. Negative staining electron micro-graphs of 
lipid DMPC vesicles (a), DMPC vesicles containing 
MHC class II-associated invariant chain transmembrane 
domain” (b) and DMPC vesicles containing SCoV E 
protein transmembrane domain (c). All vesicles were 
extruded though a 200 nm filter. 


SCoV E protein. Figure 3(d) presents the ATR-FTIR 
spectrum of the lipid CH, and CH, stretching 
modes of the aforementioned vesicles. The 


measured dichroism for the lipid CH, asymmetric 
stretching mode (centered at 2920cm'') is 
1.30 + 0.13, yielding an order parameter of 
S =0.51 + 0.11. Interestingly, the order parameter 
of the protein (S = 0.82 + 0.17) is much higher 
than that of its ordering environment, the lipid 
bilayer. The reduced order of the lipid acyl chains 
is an inherent property of the system caused by 
SCoV E protein and cannot be attributed to incom- 
plete sample deposition. The reason is that the 
exceptionally high order parameter obtained for 
the protein would not be possible under such 
circumstances. Moreover, the mosaicity of the 
membranes was measured to be below 0.01°, as 
determined by X-ray reflectivity rocking scans (not 
shown). Thus, incorporation of SCoV E protein 
transmembrane domain in the lipid vesicle 
deforms the acyl chains of the lipid bilayer. 

The global change in the structure and mor- 
phology of the lipid bilayer upon incorporating 
SCoV E protein transmembrane domain was 
examined using negative staining electron 
microscopy. While vesicles without protein (Figure 
5(a)) exhibited normal globular structure, SCoV E 
protein-containing vesicles were markedly dif- 
ferent (Figure 5(c)), exhibiting extensive tubulation 
and deformation. Vesicles containing a control 
transmembrane protein, the MHC class II-asso- 
ciated invariant chain transmembrane domain,’’ 
were similar to vesicles without protein, and did 
not exhibit any tubulation (Figure 5(b)). 
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Figure 6. Results of the global searching molecular 
dynamics simulations of the two helices forming the 
hairpin of SCoV E protein. Each structure is positioned 
according to its final rotational angle of helix I and helix 
II, and colored according to energy of the bundle as indi- 
cated by the color legend (values in kcal/mol). The gray 
lines represent a vector connecting the starting position 
for each structure and the final one. The two single 
clusters found are encircled in the Figure. 
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In summery, FTIR spectroscopy has shown that 
the incorporation of SCoV E protein in the lipid 
bilayer results in a molecular disordering of the 
lipid acyl chains. This leads to the tubulation of 
the vesicles as we have shown using electron 
microscopy. Thus, the molecular defects in the 
lipid molecule are translated into a distortion of 
the entire bilayer structure. 


A structural model of SCoV E protein 
transmembrane domain 


The evidence obtained so far points conclusively 
to the fact that the SCoV E protein transmembrane 
domain forms an unusually short transmembrane 
helical hairpin: (i) The secondary structure of the 
transmembrane domain is highly helical. (ii) More 
than 80% of the protein is embedded in the lipid 
bilayer. (iii) The helical element of the protein is 
tilted from the membrane normal by no more 
than 6.3°. (iv) Phe23 residing in the middle of the 
hydrophobic stretch of the protein is located in the 
lipid head-group region. (v) The protein contains 
a center of symmetry upon which it can be 
inverted. 

What might the hairpin structure of SCoV E pro- 
tein look like? The only deviation from hairpin 
symmetry is at Asn15 (Figure 2(b)). The groups of 
Engelman & DeGrado have recently shown that 
asparagine residues in transmembrane helices 
promote strong homo-dimerization due _ to 
H-bonding.””' In this instance there is only a 
single asparagine in the two helices and as such 
the hydrogen bonds that it might form are to a 
backbone amide group, thereby driving the two 
helices to close apposition. Alternatively, an 
asparagine residue from a different protein might 
promote dimerization, forming a pseudo four- 
helix bundle. Another interaction that might take 
place between the helices is through a salt bridge 
between Glu8 and Arg38. The presence of 
opposing charged resides at both sites is highly 
conserved in the coronaviridea (Figure 1), thereby 
substantiating the possibility of an interaction 
between them. 

In order to derive a model for the protein based 
on its deduced topology and helix tilt angles, 
global searching molecular dynamics simulations 
were employed.* Multiple starting structures 
were generated in which both helices of SCoV E 
protein (Glu7-Leu21 and Val25-Arg38) formed a 
hetero-dimer. The sizes of the helices reflected the 
high helical content of the protein. Each structure 
was distinguished by the rotational angle of the 
two helices with respect to one another, in incre- 
ments of 20°, leading to a total of 18 x 18 = 324 
structures. The tilt angle between the helices was 
set to 0°, according to the high degree of orien- 
tation found experimentally. The stability of each 
structure was then determined by a molecular 
dynamics and energy minimization trajectory. The 
resulting structures were then compared in order 
to find regions of the configuration space to which 


Figure 7. Ribbon diagram of the structural model 
derived for the transmembrane domain of SCoV E pro- 
tein. The side chains of Glu8, Asn15, Phe20, Phe 23, 
Phe26 and Arg38 as well as the carbonyl group of 
Leu31 are shown in CPK representation. Figure drawn 
with the program Ribbons.*! 


structures from different starting configurations 
converged. 

The results of the global searching molecular 
dynamics simulations are shown in Figure 6, 
depicting two clusters to which structures have 
converged. The energy of the structures and the 
two clusters show a clear preference to lower 
rotational angles of the two helices (most likely 
due to the formation of a salt bridge, as discussed 
below). Thus the cluster found at rotational angles 
of 50° and 10° (helix I and II, respectively) was 
taken to represent a model for the interactions 
between the two helices (see Figure 6). 

After the packing between the two helices form- 
ing the hairpin was obtained, the loop connecting 
the helices formed by amino acid residues Ala22- 
Val24 was constructed. The entire hairpin structure 
containing the loop was subsequently subjected to 
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energy minimization and molecular dynamics. The 
resulting model is shown in Figure (7), and exhi- 
bits three prominent features: (i) a salt bridge 
between Glu8 and Arg38, (ii) an inter-helical 
H-bond between the amide side-chain of Asn15 
and Leu31 peptide carbonyl and (iii) three phenyl]- 
alanine residues located at the tip of the hairpin 
pointing outwards forming an aromatic belt. 
Finally, based on the evidence that IBV E protein 
is palmitoylated,” it is possible to speculate that 
SCoV E protein might be as well. In contrast, 
experiments attempting to label TGEV E protein 
with palmitic acid have failed.** If present in SCoV 
E protein, palmitoylation will take place on one 
or more of the three juxtamembranous cysteine 
residues. It is difficult to estimate the effects, if 
any of palmitoylation on the SCoV E protein. 
One might speculate that the intercalation of the 
palmitic acid in the cytoplasmic membrane leaflet 
might further affect the local membrane structure. 


Relationship to in vivo studies 


In any in vitro study the relevance of the results 
to those found under native conditions is 
obviously a concern. However, the following facts 
negate any suspicions of artifactuality: 

Comparison to the full length protein. The protein 
under study encompasses all the hydrophobic 
domain of the SCoV E protein. The only thing 
known about the endogenous protein is that it 
traverses the lipid bilayer twice.’° That is exactly 
what was found herein. Thus it is highly unlikely 
that the full length protein and its membrane 
domain both traverse the lipid bilayer twice, but 
do so in different ways. 

Structural similarity to that found in native conditions. 
The endogenous function of SCoV E protein is in 
viral budding by way of causing membrane 
tubulation.*-”° Moreover, it is the only protein in the 
virus that is required for this process.”” Thus, the 
finding that the membranous segment of SCoV E 
protein (analyzed in vitro) causes membrane tubula- 
tion in an indistinguishable manner to the wild-type 
protein, in vivo, is proof that the membranous seg- 
ment is necessary and sufficient for SCoV E protein 
role in membrane tubulation. 

Finally, the membranous segment of SCoV E pro- 
tein was not studied in detergent micelles (as 
is customary), but rather in hydrated lipid bilayers, a 
system that is closest to the native conditions as poss- 
ible. It was shown that the protein has a defined 
helical secondary structure and that it is highly 
embedded in the lipid bilayer. It is difficult to think 
of an example in which a small membrane protein 
was reconstituted in lipid vesicles, shown that it was 
helical as well as highly embedded in the lipid bilayer 
but is nevertheless in a non-native conformation. 


Functional implications of SCoV E 
protein structure 


With a structural model of SCoV E protein at 


hand, is it possible to rationalize the effects of the 
protein upon the lipid bilayer? The hydrophobic 
stretch of both hairpin helices is only 11 residues 
long, which is much shorter than the average of 
21 residues.'* At the apex of the hairpin three 
amino acid residues have their amide carbonyl 
and amine groups unpaired, surrounded by an 
aromatic belt. At the other end of the hairpin a 
cluster of charged residues is located. The distance 
between these two polar elements is much shorter 
than the hydrophobic thickness of the lipid bilayer. 
The polarity of both elements is not equal however, 
whereby the charged residues at the opening of the 
hairpin predominate. Thus, one can speculate that 
the hairpin opening is situated at the lipid head- 
group region of one leaflet. Therefore, the hairpin 
apex is situated in the hydrophobic region of the 
bilayer, with its aromatic belt (known for its affin- 
ity to the polar head group) and unpaired amide 
groups poised to attract the lipid head-group region. 
It is this attraction that exerts a deformation force on 
the bilayer, causing an increase in curvature. 

In comparison to other coronaviridea E proteins, 
only the length of the hydrophobic segment lined 
by opposing charged residues at the either end, is 
conserved (Figure 1). Thus, E proteins from dif- 
ferent coronavirus subfamilies do not exhibit any 
sequence similarity to those from other sub- 
families. It would seem therefore, that members of 
the coronaviridea have “found” at least four inde- 
pendent sequence motifs to form a functional E 
protein. The effects of the unique features of SCoV 
E protein—the aromatic belt, the inter-helical 
H-bond of Asn15 and the inversion symmetry on 
the overall virus life cycle and resulting patho- 
genicity, remain to be seen. 


Materials and Methods 


Viral strains 


Abbreviations for viral strains listed in Figure 1, are as 
follows, with the accession numbers in parentheses: 
AIBV—Avian infectious bronchitis virus Beaudette US 
(CAC39303); CEC—Canine enteric coronavirus strain 
INSAVC-1 (P36696); + PTGE—Porcine transmissible 
gastroenteritis coronavirus strain FS772/70 (P2265); 
PRC—Porcine respiratory coronavirus strain RM4 
(P24415); HCV1l—Human coronavirus strain 229E 
(P19741); Equine coronavirus strain NC99 (AAQ67208); 
BCV1—Bovine coronavirus strain Mebus (AAA42914); 
BCV2—Bovine coronavirus strain F15 (P15775); 
HECV—Human_ enteric coronavirus strain 4408 
(AAQ67200); HCV2—Human coronavirus strain OC43 
(AARO1017); MHV1—Murine hepatitis virus strain A59 
(NP_068673); RSDCV—Rat sialodacryoadenitis corona- 
virus strain SDAV-681 (AAF97741); MHV2—Murine 
hepatitis virus strain S (P29076); SCoV—SARS corona- 
virus strain TOR2 (NP_828854). 


Protein synthesis, purification and reconstitution 


Peptides encompassing the predicted transmembrane 
domain of SARS coronavirus E protein (Figure 1, residues 
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7 to 38) were synthesized by standard solid-phase N-(9- 
fluorenyl) methoxycarbonyl chemistry. The peptides 
were cleaved from the resin with trifluoroacetic acid 
(TFA, Aldrich) and lyophilized. Two different synthetic 
peptides were made: an unlabeled peptide and one con- 
taining a para-iodo phenylalanine at position 23 of the 
sequence, as shown in Figure 1. Peptide purity was con- 
firmed by mass spectrometry. The lyophilized peptide 
was dissolved in TFA (final concentration ~5 mg/ml), 
and immediately injected on a Jupiter 54 C4 300 A col- 
umn (Phenomenex, Chesire, UK), equilibrated with 80% 
HO, 8% (v/v) acetonitrile (J.T. Baker and 12% (v/v) 2- 
propanol (J.T. Baker Peptide elution was achieved with 
linear gradient to a final solvent composition of 40% 
(v/v) acetonitrile, and 60% (v/v) 2-propanol. All sol- 
vents contained 0.1% (v/v) TFA. After over night lyo- 
philization of pooled fractions in the presence of 10 mM 
HCI (to remove traces of TFA adducts), ca 1 mg of the 
dried peptide was dissolved in a solution of 10mg 
dimyristoylphosphocholine (DMPC, Avanti) in 900 jl of 
1,1,1,3,3,3-hexafluoro-2-propanol (HFIP, Merck). The 
solvent was evaporated over night under reduced 
pressure, by means of a rotary evaporator. One millilitre 
of H,O was added to the dried product and the solution 
was mixed for 20 minutes at 30 °C. 


FTIR spectroscopy 


Data were recorded on a Nicolet Magna-560 infrared 
spectrometer (Nicolet Instrument Corporation, USA) 
purged with dry air and equipped with an MCTA 
detector, cooled with liquid nitrogen. A total of 1000 
interferograms were collected at a resolution of 4cm™'. 
Attenuated total reflection (ATR) spectra were measured 
with a 25 reflections ATR accessory from Grasbey Specac 
(Kent, UK) and a wire grid polarizer (0.25 ~M, Graseby 
Specac). Three hundred microlitres of sample (ca 1 mg/ml 
of peptide and 10 mg/ml of lipid) were deposited onto 
a Ge trapezoidal internal reflection element 
(50 x 2 x 10 mm) flowed by removal of bulk solvent. For 
the purpose of solvent exchange, air was followed over 
the sample that was bubbled through *H,O or Hp. 
O. Transmission FTIR spectra were collected by deposit- 
ing 50 wl of sample (ca 1 mg/ml of peptide and 10 mg/ 
ml of lipid) on a CaF, window. 

The dichroic ratios of the Amide I band were calcu- 
lated by integrating between 1670 cm™' and 1645 cm™'. 
All integrations were performed by using a straight 
baseline that contains points immediately before and 
after the band. Reported standard deviations values 
represent a minimum of three data sets measured. 


Order parameter calculation 


Order parameters for the transmembrane helices were 
calculated as described.’’ In brief, the order parameter S 
is defined as follows: 


_ 3(cos” 6) — 1 
~ 2 


whereby 0 represents the angle between the helix direc- 
tor and the z axis which is coincident with the membrane 
normal. The order parameter can be determined experi- 
mentally from the ATR dichroic ratio RATS by the follow- 
ing equation: 


S (1) 


S 


_ EZ — RATRET+E2 /3cos?a—1\" 
E2 — RATRE? — 2F2 2 
x y z 


whereby the electric field components of the evanescent 
wave for a Ge internal reflection element, E; = 1.398, 
E, = 1.516 and E, = 1.625 are given by Harrick,* and a 
is the angle between the helix director and the transition 
dipole moment of the vibrational mode: 39° for the 
amide I mode and 27° for the amide A mode.”° Lipid 
order parameters are obtained from the lipid methylene 
symmetric (2852cm™') and asymmetric (2924 cm’) 
stretching modes using the same equation by setting 
a = 90°. These equations are based on the reasonable 
assumption that the thickness of the deposited film 
(>20 ym) is much larger than the penetration depth 
(ca 1 wm) of the evanescent wave.”° 


X-ray scattering measurements 


For X-ray reflectivity measurements, samples of 
varied molar peptide-to lipid ratio (P/L= 
0, 0.002, 0.01, 0.05,0.1,0.13) were deposited on silicon 
surfaces by spreading from HFIP/chloroform solvent 
(0.1ml of 20mg/ml solution on wafers cut to 
15 X 25 mm/’). The orientational distribution (mosaicity) 
of the highly oriented membranes was measured to 
be below 0.01°, as determined by a rocking scan. High 
resolution reflectivity scans up to vertical momentum 
transfer gz ~ 1 A! were taken at the D4 bending magnet 
beamline of the DORIS storage ring at HASYLAB/DESY, 
at a photon energy of 11 keV, set by a Si(111) mono- 
chromator. The curves were measured with a fast scintil- 
lation counter (Cyberstar, Oxford instruments), using 
motorized collimating slits on the incident and exit 
arms. The curves were corrected for ring current, sample 
illumination, and diffuse background (offset-scan). 

General aspects of reflectivity experiment and analysis 
are discussed in Ref. 18 The samples were kept in an 
environmental chamber at controlled relative humidity 
of 98% and a temperature of T = 45 °C to assure that the 
bilayers were in the fluid L, state. The electron density 
profile was obtained by the Fourier synthesis method 
from the integrated peak intensities, using a Lorentz cor- 
rection factor 1/qz, and phases (+, —, —, —, —). The inte- 
grated intensity of the first Bragg peak was normalized 
to 1 for the iodinated SCoV E protein and for the Influ- 
enza A M2 curve, while the non-iodinated SCoV E 
protein curve was properly scaled to the iodinated 
curve corresponding to the respective scattering signals. 
Full q,-range fits to determine the density profiles on an 
absolute scale and analysis of all samples at different 
R.H.% is in progress (Z.K. et al., to be published). 


Electron microscopy 


Liposomes containing the transmembrane domain of 
the SCoV E protein (residues 7—38) and the MHC class 
Il-associated invariant chain (residues 29-60)'° where 
prepared as described. Liposomes composed of pure 
DMPC were prepared in an identical manner, with no 
addition of peptide. Unilamellar vesicles where obtained 
by two cycles of freeze-thaw sonication (40 seconds) 
using an inverted probe sonicator (Vibra cell model 
vcl30, Sonics and Materials Inc., CT). The liposomes 
where passed twice through 200nm_ polycarbonate 
filter using an Avanti Mini-extruder (Avanti Polar Lipids 
Inc., AL). 

Negative staining for transmission electron micro- 
scopy was done by mixing the liposome suspension at 
1:1 ratio with 2% (v/v) phosphotungstic acid (Electron 
Microscopy Sciences, Hatfield, PA) and incubating for 
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one minute at room temperature. The mixed sample was 
absorbed onto carbon-coated grids that had been ren- 
dered hydrophilic by glow discharge in air. Excess sol- 
ution was removed by filter paper absorbance. After 
drying, the samples were imaged with a TECNAT 12 
electron microscope (Philips) operating at 100 kV. 


Molecular modeling 


Molecular modeling of the transmembrane domain of 
SCoV E protein was undertaken using the CHI suite of 
macros™ for CNS 1.1.” Briefly, bundles of helices were 
generated in which the angle of the helices about their 
director was rotated independently in increments of 20°. 
The crossing angle of the helices was set close to that 
found experimentally: 0°. Each structure was subjected 
to the following simulated annealing and energy mini- 
mization protocol as follows: 600 steps of Powell energy 
minimization followed by torsion angle molecular 
dynamics at 600 K for 5 ps and at 300 K for 10 ps. Finally, 
an additional 250 steps of Powell energy minimization 
were undertaken. Throughout the simulations, harmonic 
restraints were employed to maintain helical geometry 
and distance restraints were employed to maintain the 
helix centers within 10.5 A one another.” 

The resulting 324 structures were compared, whereby 
clusters were defined as those that contained at least 
seven structures and a C* RMSD <1.0A. The clusters 
were than averaged and subjected to the same molecular 
dynamics/energy minimization protocol as above. 
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